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ABSTRACT 

We investigate the performance of spherical wavelets in discriminating between stan- 
dard inflationary models (Gaussian) and non-Gaussian models. For the later we con- 
sider small perturbations of the Gaussian model in which an artificially specified skew- 
ness or kurtosis is introduced through the Edgeworth expansion. By combining all the 
information present in all the wavelet scales with the Fisher discriminant, we find that 
the spherical Mexican Hat wavelets are clearly superior to the spherical Haar wavelets. 
The former can detect levels of the skewness and kurtosis of « 1% for 33' resolution, 
an order of magnitude smaller than the later. Also, as expected, both wavelets are bet- 
ter for discriminating between the models than the direct consideration of moments 
of the temperature maps. The introduction of instrumental white noise in the maps, 
S/N — 1, does not change the main results of this paper. 



1 INTRODUCTION 



Most of the analyses of Cosmic Microwave Background 
(CMB) data focus on the measurement of the power spec- 
trum of temperature fluctuations. Information on this sec- 
ond order moment is crucial to determine the fundamental 
parameters of the cosmological model corresponding to our 
universe. However, this determination relies on the Gaussian 
hypothesis for the temperature distribution. Establishing 
the statistical character of the CMB fluctuations will provide 
crucial evidence about the physical origin of the primordial 
density fluctuations in the early universe. Simple inflation- 
ary models predict a Gaussian, homogeneous and isotropic 
random field for the temperature fluctuations. On the con- 
trary, non-standard inflation and cosmic defects generically 
predict non-Gaussian random fields. Recent CMB observa- 
tions by Boomerang, DASI and MAXIMA-1 (Netterfield et 
al. 2001, Pryke et al. 2001, Stompor et al. 2001) have es- 
tablished for the first time the presence of multiple acoustic 
peaks in the CMB power spectrum. As a consequence cos- 
mic defects cannot be the dominant source of density per- 
turbations in the universe. Even if they are present as a 
sub-dominant component confirmation of its existence will 
be best made by appropriate techniques searching for non- 
Gaussian features in the CMB maps. 

Since a random field can departure from a Gaussian 
one in many different ways there is not a unique way to 
detect and characterise deviations from Gaussianity. Thus, 
depending on the kind of features one is looking for some 
specific methods will prove to be more efficient than others. 
Efficient methods are able to extract relevant information 
on the non-Gaussian nature of the data which is otherwise 
hidden in the temperature fluctuation maps. A large number 
of methods have been already proposed to search for non- 



Gaussianity in CMB maps. The methods can be grouped by 
the spaces (real, Fourier,...) in which they act. In real space, 
standard quantities used are the cumulants which contain 
information on the 1-pdf only. Information on the n-pdf can 
be obtained through the Edgeworth expansion (Contaldi et 
al. 2000) or alternative expansions with a proper normal- 
ization (Rocha et al. 2000). Other quantities focus on topo- 
logical and geometric statistics, e.g. Minkowski functionals 
implemented on the sphere (Schmalzing and Gorski 1998); 
statistics of excursion sets, e.g. characteristics of peaks (Bar- 
reiro, Martmez-Gonzalez and Sanz 2001), extrema corre- 
lation function (Heavens and Gupta 2001); also geometri- 
cal characteristics of polarisation have already been inves- 
tigated (Naselsky and Novikov 1998). Multifractal analysis 
and roughness have been applied to the COBE-DMR data 
(Diego et al. 1998, Mollerach et al. 1999). In Fourier space, 
the bispectrum has been applied in several occasions to anal- 
yse the COBE-DMR data (see e.g. Ferreira et al. 1998) as 
well as an extension to include possible correlations among 
multipoles (Magueijo 2000). An alternative approach is to 
work in eigen space, extracting the eigenmodes from a prin- 
cipal component analysis. This approach has been taken by 
Bromley and Tegmark (1999) for the COBE-DMR data and 
by Wu et al. (2001) for the MAXIMA-1 data. In spite of all 
this effort there is not any strong evidence of deviations from 
Gaussianity in the CMB up to date (see however Magueijo 
2000 for a possible deviation). More definitive conclusions 
about the statistical distribution of the CMB fluctuations 
are expected from data analyses of present and future sen- 
sitive experiments at arcmin resolution. 

In this work we concentrate on wavelet analyses. As it 
is often pointed out, wavelets are a very useful tool for data 
analysis due to its space-frequency localisation. It has been 
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already demonstrated in many applications in a wide vari- 
ety of scientific fields. In particular in relation to the CMB 
the COBE-DMR data has been studied with several wavelet 
bases acting on the faces of the quad-cube COBE pixelisa- 
tion (Pando et al 1998, Mukherjee et al. 2000, Aghanim et 
al. 2001). More appropriate analyses should involve the use 
of spherical wavelets as in Tenorio et al. (1999). More re- 
cently Barreiro et al. (2000) and Cayon et al. (2001) have 
convolved the COBE-DMR data with spherical wavelets in 
the HEALPiX pixelisation (Gorski, Hivon & Wandelt 1999) 
to test the Gaussianity of these data. Those works used the 
Spherical Haar Wavelet (SH W) and the Spherical Mexican 
Hat Wavelet (SMHW), respectively. 

It is our aim in this work to confront the performance 
of these two spherical wavelet bases proposed for discrimi- 
nating between standard inflationary (Gaussian) models and 
non-Gaussian models which contain artificially specified mo- 
ments (skewness or kurtosis) in the temperature distribu- 
tion. Physically motivated non-Gaussian features can en- 
ter in the CMB maps in many ways. Cosmic defects can 
produce linear discontinuities (cosmic strings, Kaiser and 
Stebbins 1986), hot spots (global monopoles, Coulson et al. 
1994) or cold and hot spots (cosmic textures, Turok and 
Spergel 1992). Non-standard inflationary models, e.g. with 
several interacting scalar fields, are expected to produce a 
qualitatively different non-Gaussianity. In particular, mod- 
els with an extra quadratic term in the potential (Linde and 
Mukhanov 1997) generate a clear signal in the third mo- 
ment (Verde et al. 2000, Komatsu and Spergel 2001). In any 
case, it is very difficult to imagine a non-Gaussian primor- 
dial model producing no significant amount of neither of the 
two low order moments. 

The paper is structured as follows. In Section 2 we 
introduce the Spherical Mexican Hat Wavelets (SMHW). 
Section 3 summarizes the main properties of the Spherical 
Haar Wavelets (SHW) and the procedure to calculate their 
coefficients. All-sky simulated non-Gaussian CMB maps at 
arcmin resolution, with a given power spectrum and artifi- 
cially specified skewness or kurtosis, are generated in section 
4. In section 5 we present optimal statistics based on the 
wavelet coefficients to get a maximum discriminating power 
between the Gaussian and non-Gaussian temperature maps. 
The main results are given in section 6 and we summarize 
the main conclusions of the paper in section 7. 



2 THE SPHERICAL MEXICAN HAT 
WAVELETS 

Future CMB missions will provide temperature data cover- 
ing all or almost all the sphere at arcmin resolution. It is 
thus necessary to have convenient pixelisation of the sphere 
which allows efficient analyses of the data. Wavelets defined 
on the plane have been widely used in astrophysical applica- 
tions during the last years. In particular, the Mexican Hat 
wavelet family has been successfully used to extract point 
sources from CMB maps (Cayon et al. 2000, Vielva et al. 
2001). However, applications of spherical wavelets have been 
very scarce and limitted to a few families of wavelets. Below 
we describe a procedure to extend the Mexican Hat wavelets 
to the sphere. 



2.1 The MEXHAT on R 2 

A continuous wavelet family on the plane R 2 is a set of filters 

built from a mother wavelet tp(x), ^(x;b,R) = -^ip( ) 
(we only consider isotropic wavelets). ip(x) satisfies the fol- 
lowing conditions: 

dxtp(x)—0 (compensation), (1) 
r 

G/, = (27r) / dqq~ tp (q) < oo (admissibility), (2) 



where b defines a translation and R a scale, i. e. we consider 
a 3-parameter family of filters. tl>(q) is the Fourier transform 
of ij) and we have introduced the standard normalization 

= = -14 (3) 



a) Analysis 

Let us consider a function on the plane f(x). The con- 
tinuous wavelet transform with respect to ^ is defined as 
the linear operation 



w(b,R)= dx f(x)^(x;b,R) 



dxf(x+b)ip(x/R). (4) 



w(b, R) are the wavelet coefficients dependent on 3 parame- 
ters. 

b) Synthesis 

It can be proven that for any tj) the following equality 
holds 



dRdbR-^(^-^)^(^—^-) = C^8(x-x'), (5) 



R 



R 



where S(x) is the Dirac distribution. 

A straightforward calculation based on the previous 
equation leads to the continuous reconstruction formula 



fix) = ^r [ dRdbR' 4 w(b,R)^-^-^-). 

J it 



(6) 



c) The MEXHAT wavelets 

A particular example is the MEXHAT wavelet defined 

by 



V(x;R) = V(x;Q,R) 



(2n) 1/2 R [2 { R 



•(7) 



This wavelet (introduced by Marr (1980)) is proportional 
to the 2D Laplacian of the Gaussian function. It has been 
extensively used in the literature to detect structure on a 
2D image (e.g. in astrophysics to detect point sources in a 
noisy background). 



2.2 The MEXHAT on S 2 

For CMB analyses we are interested in the extension of these 
isotropic wavelets to the sphere. Recently, Antoine & Van- 
dergheynst (1998) have followed a group theory approach to 
deal with this problem. This extension incorporates four ba- 
sic properties: a) the basic function is a compensated filter, 
b) translations, c) dilations and d) Euclidean limit for small 
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Figure 1. Stereographic projection to translate the properties 
of the Mexican Hat Wavelet from the plane to the sphere 



to detect non-Gaussianity in the CMB sky 3 
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Figure 2. Mexican Hat Wavelet as deformed on the sphere (solid 
line) from the plane (dashed line). The scale of the wavelet is 
chosen to be R = 1 rad. 



angles. They conclude that the stereographic projection on 
the sphere is the appropriate one to translate the mentioned 
properties from the plane to the sphere. Such a projection 
is defined by (x) t— > (0, <fi) 



x\ = 2 tan — cos ( 
2 



X2 — 2 tan - sin < 
2 



(8) 



where (0, cf>) represent polar coordinates on S 2 and (y = 
2 tan | , (fj) are polar coordinates in the tangent plane to the 
North pole (see Figure 1). 

Therefore, the isotropic wavelet ty(x;R) transforms to 



Vs(0;R) oc (cos-) 



2tan-;7?). 



(9) 



It can be proven that the new wavelet on S 2 incorpo- 
rates the basic properties, i. e. a) it is a compensated filter 
(J d6d(j} sin 89 s(d; R) = 0), b) translations are defined by 
translations along the sphere, i. e. rotations about the cen- 
ter of the sphere, c) the dilations are defined by the stereo- 
graphic projection of dilations on the plane and d) for small 
angles one recovers the Euclidean limit. 

a) Analysis 

Let us consider a function on the sphere f(8,<f>). The 
continuous wavelet transform with respect to 9s(9',R) is 
defined as the linear operation 



w(x, R) 



dO' 



'sine' f(x + p)9 s (e';R). 



(10) 



x = 2 tan — (cos <j>, sin <f>), pi =2 tan "^"( cos </>', sin </>'), 

/(x) = /(#,</>), (11) 

w(6,(j>;R) = w(x,R) are the wavelet coefficients dependent 
on 3 parameters. 

b) Synthesis 

A straightforward calculation based on the equation (5) 
leads, after stereographic projection, to the continuous re- 
construction formula: 



/(0,0 = /(S) = JL J dB'd^' sin 6'^w(x + p,R) 



9 s (e';R) 



(12) 



where w(x, R) = w(8, <j>; R). 



c) The MEXHAT wavelets 

A particular example is the MEXHAT wavelet defined 
by (see Figure 2) 



<f(8;R) = 



,4 1/2 



N(R) = (l + ^- + ^) ,y = 2tan^. 



(14) 



We remark that the normalization constant has been 
chosen such that f d6 d(f> sin 0V 2 (8; R) = 1. This is the 
wavelet we are going to use in this paper to analize non- 
Gaussianity associated to different models. 

We comment that the stereographic projection of the 
MEXHAT wavelet has been recently used to analize maps of 
the cosmic microwave background radiation (CMB) (Cayon 
et al. 2001). 



3 SPHERICAL HAAR WAVELETS 

SHW were introduced by Sweldens (1995) as a generaliza- 
tion of planar Haar wavelets to the pixelised sphere. They 
are orthogonal and adapted to a given pixelisation of the sky 
which must be hierarchical, contrary to the SMHW which 
are non-orthogonal and redundant. However they are not ob- 
tained from dilations and translations of a mother wavelet, 
contrary to planar Haar wavelets and SMHW. As for the 
planar Haar wavelets, they possess a good space-frequency 
localisation. However, their frequency localisation is not as 
good as that of the SMHW. Two applications of SHW to 
the analysis of CMB maps have already been performed. 
Tenorio et al. (1999) apply them to simulated CMB skies 
on the QuadCube pixelisation. They study the CMB spa- 
tial structure by defining a position-dependent measure of 
power. Also they show their efficiency in denoising and com- 
pressing CMB data. Barreiro et al. (2000) tested the Gaus- 
sianity of the COBE-DMR data on the HEALPix pixelisa- 
tion. One of the advantages of HEALPix over QuadCube is 
that there is no need to correct for the pixel area. 
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Figure 3. Hierarchical structure of Wavelet coefficients for the 
Spherical Haar Wavelet. 



Since detailed description of the SHW transform has al- 
ready been given in the previous papers, we here describe the 
main features of the wavelet decomposition. The SHW de- 
composition is based on one scaling <j>j t k and three wavelet 
functions ip m ,j,k at each resolution level j and position on 
the grid k. For HEALPix the resolution is given in terms of 
the number of divisions in which each side of the basic 12 
pixels is divided, N si d e — 2 3 . Thus, for level j the total 
number of pixels with area fij is given by rij — 12 x 4P~ . 
Each pixel k at resolution j, Sj t h is divided into four pix- 
els Sj+i t k Sj+i,k 3 at resolution j + 1. For computational 
reasons the maximum resolution we will consider in our sim- 
ulations is J = 9 which corresponds to N s ide ~ 256. The 
scaling and wavelet functions are simply given by 



ipi,j,k = 

1p2,j,k = 
Tp3,j,k = 



i ifx e s 3]k 

otherwise 



<f>i+l,k 


+ 












-1 


0i+l,fco 


+ 


Pj+lM - 


<t>j+-LM - 0j+l,fc3 








-i 




+ 


t>i+l,k 3 - 


4>j+iM — <j>j+i,k 2 



(15) 
(16) 
(17) 
(18) 



where ko,ki,k2, k$ are the four pixels at resolution level j + 1 
in which the pixel k at level j is divided. Please note that 
the three wavelet functions so defined differ from the ones 
used by Tenorio et al (1999) and Barreiro et al. (2000). We 
choose those expressions by similarity with the diagonal, 
vertical and horizontal details defined on the plane. The re- 
construction of the temperature field is obtained by 

"JO" 1 J-IUj-1 

{Xi)= ^ ^o^MJ^l+Xil] ^ 7m,j,i 

(=0 m j=j 1=0 

tpm,j,l(xi) , (19) 



AT 



where Aj 0i fc and 7m,j,fc are the approximation and detail co- 
efficients respectively. The level index j goes from the finest 
resolution J to the coarsest one considered jo. 

The wavelet coefficients at level j can be obtained from 
the four corresponding approximation coefficients at level 
j + 1, Aj+i i fc i as follows (see figure ^.): 



3 



(20) 





i — 


ll 












= N- 




+ ^j+l,k 2 


— Aj + l,*;! 


- A J+ i,fc 3 ) 


(21) 


72,j,fc 


= 


n(Aj + i,fc 


+ Aj + l.fc! 


— Aj + l,fe 2 


— \ j+ i } k 3 ) 


(22) 


73,j,k 


= Mi- 


n(^i+i,fc 


+ ^j+l,k 3 


— Aj + l^fej 


— \j+i,k 2 ) 


• (23) 



The generation of coefficients start with the original map, 
finest resolution j = J, for which the coefficients Aj.fc are 
identified with the temperature fluctuation at pixel k. 

Finally, from the definition of the SHW it is easily seen 
that this wavelet is not rotationally invariant, contrary to 
the SMHW. 



4 NON-GAUSSIAN SIMULATIONS 

There are many ways in which physically motivated non- 
Gaussian features can enter in the CMB temperature dis- 
tribution. However, up to now there is no evidence of their 
existence, being all experimental data consistent with Gaus- 
sianity (Kogut et al. 1996, Barreiro et al. 2000, Aghanim et 
al. 2000, Cayon et al. 2001, Wu et al. 2001; see however 
Magueijo 2000 for a possible positive signal in the COBE- 
DMR data, although that detection has not been confirmed 
by any of the other analyses). If departures from Gaussian- 
ity of cosmological origin really exist they will more likely 
be small and all-sky, sensitive, arcminute resolution experi- 
ments will be needed for their detection. 

Here the spherical wavelets will be tested against non- 
Gaussian simulations of artificially specified moments that 
will be assumed to be small. In this case a useful way to 
construct non-Gaussian distributions is by perturbing the 
Gaussian one through a sum of moments, the Edgeworth 
expansion. For simplicity we will consider the two lowest 
cumulants to characterise the deviations from normality: 
skewness and kurtosis. As discussed in the introduction al- 
ternative models to standard inflation, e.g. cosmic defects 
as a subdominant source of density perturbations or non- 
standard inflation, can produce significant levels of at least 
one of the two moments. 



4.1 Edgeworth expansion 

For small deviations from Gaussianity, there is a wide class 
of distributions that can be given in terms of a Gaussian 
distribution times an infinite sum of its cumulants. This is 
the well known Edgeworth expansion. The problem with this 
expansion is that setting all cumulants to zero except one 
does not guarantee the positive definiteness and normaliza- 
tion that a distribution has to satisfy. However, for small 
deviations from normality the resulting function is always 
positive at least up to many sigmas in the tail of the distri- 
bution and the normalization factor required for the function 
to become a well defined distribution is very small and does 
not appreciably disturb the non-zero moments (i.e. skewness 
or kurtosis) introduced in the first place. 

The Edgeworth expansion can be obtained from the 
characteristic function <j>(t) by considering the linear terms 
in the cumulants and performing the inverse Fourier trans- 
form to recover the density function f(x): 
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f {x) = G( X ){l + J2-^H n (^)+0(k n k n ,)} , (24) 

n — 3 

where H n is the Hermite polynomial. Considering the per- 
turbations corresponding to the skewness and kurtosis and 
keeping only the first terms in the corresponding Hermite 
polynomials, we have 

fsix) = ^(l+ f(*(* 2 -3))), (25) 

x 2 

f K ( x ) = ^(l + ^(^_6x 2 + 3)) , (26) 

where S, K denote skewness and kurtosis, respectively. We 
will use these equations to generate our artificially specified 
non-Gaussian distributions. Since the resulting distribution 
is not well defined even for the case of small skewness and 
kurtosis we set the function to zero when it becomes nega- 
tive and we also normalize it appropriately. We remark that 
the zero cuts of the distribution, if present, appear far away 
in the tails of the distribution for the case of small values of 
skewness and kurtosis that we consider here. Also, as a con- 
sequence, the normalization value required is very close to 1. 
In this way we checked that the initial values of the skewness 
and kurtosis we start with in the Edgeworth expansion does 
not appreciably change after the necessary changes intro- 
duced to obtain a well defined probability density function 
(pdf). 

In order to make the simulations resemble the CMB 
data observed by a given experiment we smooth them with 
a Gaussian filter. For practical reasons we use a FWHM of 
33' which may correspond to some of the channels in all-sky 
experiments like MAP and Planck (e.g. the 30GHz channel 
of the Planck mission) . We choose to work on the HEALPix 
pixelisation of the sphere with a resolution N B id e = 256. We 
use the HEALPix package to perform the analysis of our 
simulated CMB data. However, it is not adequate to use 
that package to convolve our unfiltered independent temper- 
ature data with the Gaussian 33' FWHM beam in Fourier 
space, instead we perform the convolution in real space. Af- 
ter that, in order to make the simulations more realistic 
we normalize the CMB power spectrum Ci of both Gaus- 
sian and non-Gaussian simulations to that of a CDM flat A- 
model using the HEALPix package. As a consequence of the 
beam convolution and the introduction of correlations in the 
temperature maps the original levels of skewness and kur- 
tosis injected through the Edgeworth expansion are reduced 
(compare columns 1 and 2 in table |^). The performance of 
spherical wavelets will be tested with these simulations in 
section 5. 

4.2 Distribution of spherical wavelet coefficients 

Since wavelet coefficients represent linear transformations 
of the original data, in the case of a Gaussian distribution 
the wavelet coefficients remain Gaussian distributed. This a 
very nice property of wavelets and all we have to do to test 
Gaussianity in wavelet space is to look from deviations from 
normality. 

However, for the case of the sphere any given pixeli- 
sation scheme will introduce biases. The specific bias in- 
troduced will depend on, for instance, whether the pixels 



are not of equal area or the distances between one pixel 
and its neighbours vary with the position on the sphere. 
This is in fact the situation for the two pixelisations already 
used to analyse all-sky CMB temperature fluctuations. For 
the COBE-DMR experiment the pixelisation used was the 
Quad-Cube and in this projection of the cube on the sphere 
equal-area pixels on the sides of the cube appear with differ- 
ent area when projected on the sphere. For present satellite 
experiments like MAP and Planck the HEALPix pixelisa- 
tion is now widely used. While this pixelisation possesses 
very nice properties, such as equal area iso-latitude pixels, 
however the distances between one pixel and its neighbours 
vary with latitude. Pixels near the equator tend to be more 
uniformly distributed than those near the poles. As we will 
compute in next section, this property produces a bias in 
the kurtosis of the wavelet coefficients for the case of the 
SHW (see table 1, Gaussian case which corresponds to a 
null injected value for the kurtosis). For the Gaussian and 
non-Gaussian simulations which will be performed in next 
section we will compute the first cumulants of the coefficients 
of the two spherical wavelets considered in this paper for the 
HEALPix scheme. For the SHW the coefficients correspond 
to three different details: diagonal, vertical and horizontal. 
Since those details are directly obtained from linear opera- 
tions of the four neighbour pixels (as we saw in the previous 
section) and pixels are not equally separated all over the 
sphere, correlations present in the temperature fluctuations 
make the wavelet coefficients to be biased. This bias pro- 
duces a peaked distribution with respect to a Gaussian and 
therefore a positive kurtosis in the three details of the SHW 
coefficients even for temperature realizations derived from 
normal distributions (as can be seen from table 1, the mean 
value of the kurtosis for the finest resolution of the Gaussian 
model is displaced about 10<r from zero). 

In the case of the SMHW we only have a type of coeffi- 
cients for each scale. Since this is a continuous, rotationally 
invariant wavelet -and thus not adapted to the pixelisation- 
no bias is produced in this case. 



5 DISCRIMINATING POWER 

The discriminating power of the spherical wavelets will 
be tested using Gaussian and non-Gaussian simulations 
with different amounts of either skewness or kurtosis in- 
troduced using the Edgeworth expansion, and normalized 
to a power spectrum Ci consistent with observations (as 
discussed above). Since the skewness and kurtosis are in- 
troduced at the highest resolution through the Edgeworth 
expansion (as described above), we expect to detect them 
with the skewness and kurtosis of the spherical wavelet co- 
efficients also at the highest resolutions. Thus we will con- 
sider for the analysis the first five resolution scales starting 
from the finest one. The scales go as powers of 2 for the 
SHW and for comparison we choose the same values for the 
SMHW parameter R: 1, 2, 4, 8 and 16 pixels. We can relate 
the scales of the two wavelets by looking to the scaling func- 
tions. The relation between the side, s, of the step function 
(scaling function for the Haar wavelet) and the dispersion 
R of the Gaussian is: s = y/2nR. Then, for the finest scale 
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Table 1. Mean and standard deviation, within parenthesis, for the different wavelet scales. 
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s = 2 pixels, which corresponds to an R ~ 0.8 pixels which 
is approximately 1 pixel. 

Results obtained in Fourier space are equivalent to those 
obtained in real space if the functions considered are band- 
width limited (with the bandwidth included in the one cov- 
ered by the pixelisation) . We have checked this for the finest 
resolution of the SMHW. The average difference between 
the SMHW coefficients computed by direct convolution in 
real space and going to Fourier space is < 1%. 

Given the 5 values of skewness or kurtosis corresponding 
to the 5 resolution scales for the SMHW and the 15 values 
for the SHW (5 scales for each of the 3 details), we would 
like to construct a test statistic which, combining all this in- 
formation, can best distinguish between the two hypotheses: 
a) Ho: the data are drawn from a Gaussian model, b) Hf. 
the data are drawn from a non-Gaussian model with either 
skewness or kurtosis. The best test statistic in the sense of 
maximum power for a given significance level is given by the 
likelihood ratio: 



t{x) 



f(x\Ho) 



(27) 



f(x\Hi) 

where f(x\Ho) and f(x\Hi) are the pdf of the data given 



hypotheses Ho and Hi, respectively. Since we do not know 
those multivariate pdfs and would be tremendously costly 
in cpu time to determine them by Monte Carlo simulations, 
we use as test statistic the simpler Fisher linear discrimi- 
nant function (Fisher 1936; see also Cowan 1998). This dis- 
criminant has been recently used by Barreiro and Hobson 
(2001) to study the discriminanting power of planar wavelets 
to detect non-Gaussianity in the CMB in small patches of 
the sky. The Fisher discriminant is a linear function of the 
data that maximizes the distance between the two pdfs, 
g(t\Ho) and g(t\Hi), such a distance defined as the ratio 
(to — ti) 2 / (Eg + Ei). Tk and E 2 ,, k — 0, 1, are the mean and 
the variance of g(t\Hk), respectively. The Fisher discrimi- 
nant is given by: 



t(x) = (fa-fii) W- 



(28) 



with W = Vo + Vi and Vk the covariance matrix and /2t the 
mean values of f(x\Hk)- In the particular case that f(x\Ho) 
and f(x\Hi) are both multidimensional Gaussians with the 
same covariance matrix, the Fisher discriminant is equiva- 
lent to the likelihood ratio. 

The mean values and covariance matrices of the skew- 
ness and kurtosis at each resolution level for the Gaussian 
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Table 2. Power of the Fisher discriminant at 1% significance level 





Injected 


True 1 


SMHW 


SHW 


Temperature 






X10~ 2 


P(%) 


P(%) 


P(%) 




0.05 


0.9(2.4) 


66.8 


1.51 


2.51 




0.10 


1.6(2.3) 


100 


7.09 


4.67 


SKEWNBSS 


0.30 


4.6(2.4) 


100 


36.12 


36.85 




0.30 2 


1.1(0.7) 


99.6 


1.80 


2.83 




0.50 


6.9(2.4) 


100 


78.46 


73.6 




0.10 


0.3(2.6) 


15.35 


3.00 


1.42 




0.30 


0.8(2.7) 


86.89 


9.00 


3.40 


KURTOSIS 


0.40 


1.1(2.7) 


98.10 


16.11 


4.90 




0.50 


1.4(2.6) 


99.90 


28.43 


3.50 




0.50 2 


0.2(0.7) 


20.84 


1.00 


0.32 



True refers to the mean value obtained in the analysed maps. The standard 
deviation is given within parenthesis. 
2 These models include the addition of noise to the maps with S/N — 1. 

and non-Gaussian models are obtained from a large number 
of simulations. In the next section we use those simulations 
to compare the power of the test p = 1 — f3 to discriminate 
against the alternative hypothesis Hi at a given significance 
level a for the two spherical wavelets, a and j3 account for 
the probability of rejecting the null hypothesis Ho when it 
is actually true (error of the first kind) and the probability 
of accepting Ho when the true hypothesis is Hi and not Ho 
(error of the second kind) , respectively. The decision to ac- 
cept or reject Ho is done by defining a critical region for the 
statistic t; if the value of t is greater than a cut value t C ut 
the hypothesis Ho is rejected. Thus, a and f3 are given by: 



m = / dtg{t\Ho) , 

'taut 

I - j ^ dtg{t\Hi) . 



(29) 
(30) 



This kind of analysis is very much along the lines of the 
one performed by Barreiro and Hobson (2001) for planar 
wavelets. From now on a value for the sensitivity of a = 1% 
will be adopted. 



6 RESULTS 

For both, Gaussian and not Gaussian models, we perform 
a 1000 simulations. As commented above, to make the sim- 
ulations more realistic each simulation is convolved with a 
Gaussian filter of 33'. In addition, its power spectrum Ci 
is normalized to that of a CDM flat A-model using the 
HEALPix package. For each of the simulations the wavelet 
coefficients for both the SMHW and the SHW are com- 
puted. The SMHW coefficients are computed by convolv- 
ing the CMB map with the SMHW given in eq. (13). We 
again use the HEALPix package to perform such convolu- 
tion in Fourier space, having previously calculated the Leg- 
endre coefficients of the SMHW at the specified resolution. 
The SHW detail coefficients are computed by performing 
the linear combinations of 4 pixels as described in section 
3. Computation time of wavelet coefficients using HEALPix 
scale as Nx side and N$f side for SMHW and SHW respec- 
tively. 

In figure 4 we show the mean values and dispersion of 
the skewness and kurtosis of the Gaussian and non-Gaussian 



models for the temperature map, and for the first 5 resolu- 
tion levels of the SHW diagonal, vertical and horizontal coef- 
ficients and SMHW coefficients. As expected the differences 
are best seen in the finer resolutions. It is clear from figure 4 
that the differences in the skewness for the two models are 
more remarkable for the SMHW than for the SHW and the 
temperature map. This is also the case for the kurtosis. As 
we pointed out in section 4.2, there is a strong bias in the 
kurtosis of the three details of the SHW coefficients due to 
the slight non-uniform distribution of pixels on the sphere 
in the HEALPix pixelisation. This kind of bias is expected 
for any pixelisation of the sphere due to the impossibility of 
having a uniform pixelisation. The specific bias introduced 
will depend on the pixelisation scheme used. On the con- 
trary, no bias is present for the SMHW coefficients due to 
its continuous nature. The Fisher discriminant t can still be 
applied to distinguish between the two models even in the 
presence of that bias in the kurtosis. As seen in the previous 
section, what enters in the linear coefficients to compute the 
statistic t is the difference between the means from the two 
models, cancelling out the bias term. In figures 5,6 we show 
the pdfs of the statistic t for three values of the skewness 
and kurtosis of the non-Gaussian models. It is clear that 
for both non-Gaussian models, with either positive skew- 
ness or kurtosis, the SMHW is able to distinguish between 
the Gaussian and non-Gaussian models much better than 
the SHW. 

In table |^ the power p of the Fisher discriminant con- 
structed from the skewness or kurtosis of the SMHW, SHW 
and temperature is given for several values of the cumulants. 
For the case of the temperature of the map the statistic is 
given directly by its cumulants. Again, the performance of 
the SMHW is superior to the SHW and the temperature in 
all cases. 

Since the SHW is affected by the non-uniform pixelisa- 
tion of the sphere, one might wonder if its failure to detect 
non-Gaussianity is a feature of the Haar wavelet in general 
or a consequence of the pixel dependent scale mixing. In 
order to answer this question we have made the same com- 
parison between Gaussian and non-Gaussian models, one 
with skewness 0.3 and the other with kurtosis 0.3, but now 
on the plane. (We have considered simulated 12 .8 x 12 .8 
maps with l'.5 pixels and a beam of 5' FWHM. The steps of 
the simulation and analysis are the same as for the sphere). 
The result is very similar to the one found on the sphere. 
Therefore, the failure of the Haar wavelet to detect non- 
Gaussianity is an intrinsic caracteristic of this wavelet and 
not a consequence of the pixel dependent scale mixing due 
to its implementation on the sphere. (Notice, however, that 
its performance can be similar to other planar wavelets for 
some specific features more adapted to its shape, e.g. cosmic 
strings, see Barreiro and Hobson 2001). The pixel dependent 
scale mixing basically induces a bias which has been taken 
into account in the analysis. 

In order to know the effect of instrumental noise (white) 
on the discriminating power of the spherical wavelets we 
have also added noise to the temperature maps with an am- 
plitude equals to the signal, S/N = 1. In this case 500 sim- 
ulations were generated. As shown in figure 7, the first reso- 
lution scale is the most affected and now the second scale is 
the most relevant for discrimination between models. In this 
figure it is also plotted the new pdfs for the Fisher discrim- 
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Figure 4. Comparison of Mexican Hat wavelet (black circle) and the Haar Wavelet details: Vertical (diamond), Diagonal (triangle) and 
Horizontal (square); for Kurtosis (left) and Skewness (right) values of 0.1, 0.3 and 0.5 (from top to bottom). Each point represents the 
number of sigmas deviated from the Gaussian model. Also plotted is the stripe for the non-Gaussianity determined from the temperature 
map (in this case only the pixel scale is meaninful, the stripe is drawn only for illustrative purposes). 



inant for 30% injected skewness and 50% injected kurtosis. 
The noise effect is shown in the narrowing of the separation 
between distributions as compared to the no-noise case. We 
see that the SMHW is still able to discriminate with a high 
power for the skewness model with a skewness value in the 



analysed map of 1.1%. For the kurtosis model, the addition 
of noise with the same amplitude than the signal reduces 
the level of kurtosis in the analysed map from 1.4% to 0.2%, 
a level too low to be detectable. 

Finally, even if future experiments like MAP and Planck 
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Figure 5. Fisher Discriminant for Skewness from Spherical Haar Wavelet (left) and Spherical Mexican Hat Wavelet (right). From top 
to bottom, the values of skewness in the original maps are 0.1, 0.3 and 0.5. The solid line is the Gaussian model, while the dashed one 
represents the non-Gaussian case. 



observe the full sky probably only the fraction outside the 
Galactic plane will be used to test non-Gaussianity. This 
problem has already been considered in previous papers 
analysing the Gaussianity of the COBE-DMR data with the 
SHW and the SMHW (Barreiro et al. 2000, Cayon et al. 
2001). As can be seen from those papers the impact on the 
two methods is similar (in both cases one looses all the co- 
efficients computed from pixels intersecting the cut). In any 



case, for future missions like MAP or Planck, the Galac- 
tic cut should be much smaller than for COBE because of 
the much better resolution and the much larger frecuency 
information, implying a smaller impact on the analysis. 
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represents the non-Gaussian case. Please note the bias in the distribution of the SHW kurtosis as discussed in the text. 



7 CONCLUSIONS 

We have compared the performance of the two spherical 
wavelet families already used to test the Gaussianity of the 
COBE-DMR CMB data: Mexican Hat (Cayon et al. 2001) 
and Haar (Barreiro et al. 2000). As testbed we use non- 
Gaussian simulations of all-sky arcminute resolution CMB 
maps, with a power spectrum C'i consistent with observa- 



tions and artificially specified amounts of skewness or kur- 
tosis. Most, if not all, physically motivated non-Gaussian 
primordial models of structure formation proposed in the 
literature show some amounts of either of these two mo- 
ments in the CMB maps. These simulated sky maps are 
pixelised using the widely used HEALPix package. As com- 
mented in section 4.2 any pixelisation scheme of the sphere 
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Figure 7. Top graphics: Comparison of Mexican Hat wavelet (black circle) and the Haar Wavelet details: vertical (diamond), diagonal 
(triangle) and horizontal (square); for skewness 0.3 (right) and kurtosis 0.5 (left) with added noise. Each point represents the number of 
sigmas deviated from the Gaussian model with noise. Also plotted is the stripe for the non-Gaussianity determined from the temperature 
map (in this case only the pixel scale is meaninful, the stripe is drawn only for illustrative purposes). Bottom graphics: Fisher discriminant 
for skewness (right) and kurtosis (left) from Spherical Haar and Mexican Hat Wavelets 



will introduce a bias because of the impossibility of a uni- 
form pixelisation. In particular, for the HEALPix scheme 
this bias shows up as a positive kurtosis in the Spherical 
Haar Wavelets (SHW) coefficients even for temperature re- 
alizations derived from normal distributions. The bias rep- 
resents a w 10a effect for the finest resolution, as can be 
seen from the first row of the kurtosis in table 1. On the 



contrary no bias is present in the case of the Spherical Mex- 
ican Hat Wavelet (SMHW) due to its continues nature, i.e. 
not adapted to the pixelisation scheme. 

The main conclusion of this paper is that the SMHW 
bases are much more efficient to discriminate between Gaus- 
sian and non-Gaussian models with either skewness or kur- 
tosis present in the CMB maps than the Spherical Haar 
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Wavelet (SHW) ones. More specifically, the SMHW is able 
to discriminate a 1.6% skewness with a power of 100% at the 
1% significance level whereas the SHW can weakly discrimi- 
nate a 6.9% skewness with a power of only 78% at the same 
significance. In the case of kurtosis, the SMHW detects a 
1.4% level with a power of 99.9% whereas for the SHW the 
power is only 28%, at the same 1% significance level. The 
failure of the Haar wavelet to detect non-Gaussianity is not 
a consequence of the pixel dependent scale mixing due to its 
implementation on the sphere but an intrinsic property of 
this wavelet (as has been demonstrated by performing a sim- 
ilar analysis on the plane). If we were to use the temperature 
map instead of the wavelet coefficients, the power would be 
always smaller than for the wavelets (only comparable to the 
SHW in the case of skewness). An interesting property of the 
SMHW is that an injected skewness/kurtosis in the temper- 
ature maps produces an amplified skewness/kurtosis in the 
SMHW coefficients and a negligible kurtosis/skewness. On 
the contrary the SHW is not able to amplify any injected 
skewness/kurtosis with neither skewness nor kurtosis of its 
coefficients. 

Finally, we have also tested the performance of the 
spherical wavelets in the more realistic case in which instru- 
mental noise (white) is present in the maps. In this case the 
highest resolution scale is the most affected, being the best 
scale for discrimination the second one. For a signal to noise 
ratio S/N = 1, and combining all the information from all 
the scales with the Fisher discriminant, the SMHW is still 
able to discriminate with a high power levels of skewness 
and kurtosis above 1%. 
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